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ABSTRACT 


This  paper  considers  a population  (composed,  for  example,  of  people, 
machines,  or  livestock)  subject  to  a randomly  occurring  defect  or  disease. 
If  there  exist  testing  procedures  capable  of  detecting  the  defect  before 
it  would  otherwise  become  known,  and  if  such  early  detection  provides  bene- 
fit, the  periodic  administration  of  such  a test  procedure  to  the  members 
of  the  population,  i.e.,  a mass  screening  program,  may  be  advisable.  This 
paper  develops  and  analyzes  two  important  cases  of  a general  model  of  a 
mass  screening  program.  In  both  cases,  the  defect  arrival  process  is 
Poisson. 

In  the  first  case,  the  reliability  of  the  test,  i.e.,  the  probability 
that  the  test  will  detect  the  disease,  is  dependent  only  on  the  type  of 
test  being  administered  and  is  constant  over  all  values  of  elapsed  time 
since  the  incidence  of  the  defect.  The  long-run  expected  disutility  per 
unit  time  is  derived  and  variations  of  the  disutility  with  regard  to  test 
reliability  and  testing  frequency  are  presented.  In  addition,  explicit 
solutions  are  provided  for  various  special  forms  of  the  disutility  func- 
tion. 


In  the  second  case,  the  test  will  detect  the  defect  if  and  only  if 
the  elapsed  time  since  incidence  exceeds  a critical  threshold  T which 
characterizes  the  test.  An  expression  for  the  long-run  expected  disutil- 
ity per  unit  time  is  derived  and  decision  rules  to  select  between  two  al- 
ternative test  types  and  costs  are  presented. 


• * ! Si 

! /I  w 

3 ? * 4 

9 » 3 § I r.  U ^4  j 

V m '•  < 8 s - I**,  j 


OPTIMAL  MASS  SCREENING  UNDER  CONSTANT  AND 
THRESHOLD  TEST  RELIABILITIES 


Mass  screening  is  the  process  of  inspecting  a large  population  for 
defects.  This  population  may  consist  of  many  subpopulations  each  with 
different  numbers  of  members  but  more  importantly  each  with  different 
levels  of  susceptibility  to  the  defects.  If  the  early  detection  of  a 
defect  provides  benefit,  it  may  be  advantageous  to  employ  a test  capable 
of  revealing  the  defect's  existence  in  its  earlier  stages.  (Throughout 
this  paper  the  words  "defect"  and  "unit"  or  "individual"  will  refer  to 
defect,  disorder  or  disease  and  to  a member  of  the  population,  respec- 
tively. ) 

Defects  may  arrive  in  a seemingly  random  fashion  such  as  many  types 
of  machine  failure,  the  incidence  of  certain  types  of  cancer,  diabetes, 
glaucoma,  heart  disease,  etc.,  or  they  may  arrive  as  the  result  of  some 
contagion  such  as  smallpox,  polio,  etc.  It  is  the  former  type  of  arrival 
process,  random  arrivals,  which  is  studied  in  this  paper. 

Now  continuous  monitoring  would  provide  the  most  immediate  detection 
of  a defect.  But  considerations  of  expense  and  practicality  will  fre- 
quently rule  out  continuous  monitoring  so  that  a schedule  of  periodic 
testing--a  screening  program--may  be  the  most  practical  means  of  achieving 
early  detection  of  the  defect.  In  general  terms,  the  question  then  be- 
comes one  of  how  best  to  trade  off  the  expense  of  testing  which  increases 
both  with  the  frequency  of  test  applications  and  with  the  cost  of  the 
type  of  test  used  against  the  benefits  to  be  achieved  from  detecting  the 
defect  in  an  earlier  ctage  of  development. 
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The  benefits  of  early  detection  depend  upon  the  application  consi- 
dered. For  example,  in  a human  population  being  screened  for  some  chronic 
disease  the  benefits  of  early  detection  might  include  an  improved  proba- 
bility of  ultimate  cure,  diminished  time  period  of  disability,  discomfort, 
and  loss  of  earnings,  and  reduced  treatment  cost.  If  the  population  being 
screened  consists  of  machines  engaged  in  some  kind  of  production,  the  bene- 
fits of  early  detection  might  include  a less  costly  ultimate  repair  and  a 
reduction  in  the  time  period  during  which  a faulty  product  is  being  un- 
knowingly produced.  If  the  population  being  screened  consists  of  machines 
held  in  readiness  to  meet  some  emergency  situation,  an  early  detection  of 
a defect  would  reduce  the  time  the  machine  was  not  serving  its  protective 
function. 

The  expense  of  testing  includes  easily  quantifiable  economic  costs 
such  as  those  of  the  labor  and  materials  needed  to  administer  the  testing. 
However,  there  can  also  be  other  important  cost  components  which  are  more 
difficult  to  quantify.  For  example,  in  the  case  of  a human  population 
subject  to  medical  screening,  the  cost  of  testing  includes  the  inconven- 
ience and  possible  discomfort  necessitated  by  the  test;  the  cost  of  false 
positives  which  entails  both  emotional  distress  and  the  need  to  do  unnec- 
essary follow-up  testing;  and  even  the  risk  of  physical  harm  to  the  testee, 
e.g.,  from  the  cumulative  effect  of  X-ray  exposure. 

In  an  earlier  paper  by  Pierskalla  and  Voelker  [1976],  a mathematical 
model  of  a mass  screening  program  over  Q subpopulations  was  developed. 

The  model  allows  for  various  testing  technologies,  with  different  relia- 
bility characteristics  and  costs,  different  arrival  rates  of  the  defect 
for  each  subpopulation,  and  budget,  personnel  and  facility  constraints. 
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The  objective  is  to  minimize  the  long-run  disutility  per  unit  time  for 
the  population.  The  decision  variables  are  the  testing  frequencies  for 
each  subpopulation  including  the  consideration  whether  a mass  screening 
program  is  justified  at  all  for  any  given  subpopulation. 

In  this  paper,  this  model  is  very  briefly  presented  and  then  used 
to  consider  two  important  special  cases  of  mass  screening  situations. 

The  first  case  is  that  of  a test  or  tests  which  have  constant  test  relia- 
bility. This  case  is  often  one  in  which  the  testing  instrument  (machine, 
vaccine,  etc.)  has  a relatively  constant  probability  of  false  negatives 
or  is  one  where  the  test  may  be  perfectly  reliable  but  the  machine  or  de- 
vice randomly  produces  defective  units  with  a relatively  constant  proba- 
bility. 

The  second  special  case  of  the  general  model  is  that  of  a test  or 
tests  which  have  a threshold  reliability.  In  this  case,  there  is  an 
elapsed  time  T i 0 after  the  arrival  of  the  defect  such  that  the  test  is 
not  able  to  detect  the  defect  up  to  time  T but  will  certainly  detect  the 
defect  once  the  defect  has  matured  beyond  time  T.  This  case  is  common  in 
practice  where  there  is  a threshold  of  the  defect  before  the  test  can 
recognize  its  existence.  This  threshold  time  T is,  of  course,  a function 
of  the  test  being  administered. 

Before  proceeding  to  the  cases  of  constant  and  threshold  test  relia- 
bilities, it  is  useful  to  briefly  discuss  some  of  the  related  literature 
on  mass  screening  models. 

Some  models  of  mass  screening  have  used  a multi-state  approach.  The 
Shwartz  and  Galliher  [ 1975]  model  was  concerned  with  mass  screening  for 
breast  cancer.  They  utilized  the  Bross-Blumenson  model  for  the  progressive 
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evolution  of  the  disease  through  a set  of  states  specific  to  breast  cancer. 
Detection  of  the  disease  carries  a certain  associated  increase  in  life  ex- 
pectancy which  is  determined  by  the  state  of  the  disease  at  the  time  of 
detection.  Moreover,  both  the  probabilities  of  detecting  the  disease  by 
means  of  a screening  test  and  of  detecting  it  in  some  other  way  (e.g.,  the 
woman  noticing  a lump)  are  a function  of  the  disease  state.  Upon  calibrating 
the  parameters  of  their  model  to  empirical  data,  Shwartz  and  Galliher  utilize 
computer  simulation  to  determine  the  payoff  due  to  particular  longitudinal 
schedules,  i.e.,  due  to  particular  sequences  of  testing  intervals  for  a 
woman  over  her  adult  lifetime. 

In  Voelker  [ 1 976] , multi-state  models  were  considered  in  which  there 
are  countably  many  defect  states  through  which  the  disorder  evolves  accord- 
ing to  a semi-Markov  process.  (Zelen  and  Weiss  [1965]  have  shown  a semi- 
Markov  process  to  provide  a good  fit  to  the  course  of  development  of  one 
disease  which  they  studied  in  detail.)  The  probability  of  a test  detecting 
an  existing  defect  is  a function  of  the  state  of  the  defect  at  the  time  the 
test  is  applied.  Certain  defect  states  are  regarded  as  "self-detecting". 

The  reward  for  detecting  a defect  through  the  use  of  a screening  test  is 
a function  of  the  state  of  the  defect  at  the  time  of  the  detection  and  how 
long  the  defect  has  been  in  that  state.  The  semi -Markov  process  of  defect 
evolution  is  independent  of  the  time  of  incidence  of  the  corresponding  de- 
fect which  is  a Poisson  process. 

The  results  include  formulas  for  the  long-run  expected  utility  per 
unit  time  of  a screening  program,  and  several  relationships  involving  the 
"steady-state"  distribution  of  undetected  defect  states  which  will  arise 
in  response  to  that  screening  program. 
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Most  works  in,  and  related  to,  the  area  of  modeling  mass  screening 
have  not  utilized  a multi-state  approach,  but  rather  a time  based  approach. 
Such  models  usually  focus  upon  either  the  detection  delay  or  the  lead  time 
characteristics  of  a screening  program.  Detection  delay  is  the  time  between 
the  incidence  of  the  defect  and  its  detection  (whether  that  detection  is 
the  result  of  a screening  test  or  of  the  defect  becoming  self-evident). 

Lead  time  is  the  time  difference  between  the  time  of  detection  via  a 
screening  test  and  the  time  detection  would  otherwise  have  occurred  had 
not  a screening  program  been  in  existence.  These  time  based  models  some- 
times do  make  reference  to  defect  states,  but  the  number  of  states  is  con- 
fined to  two,  the  preclinical  state  and  the  clinical  state.  In  the  pre- 
clinical  state,  the  defect  exists  but  has  not  been  discovered.  In  the 
clinical  state^,  the  defect  has  become  self-evident--thus  rendering  a 
screening  test  superfluous.  When  the  defect  occurs,  it  is  first  in  the 
preclinical  state  and  then  in  the  clinical  state.  Hence,  both  detection 
delay  and  lead  time  are  determined  by  the  time  of  detection  relative  to 
the  sojourn  interval  of  the  preclinical  state. 

Some  early  papers  on  a time-based  approach  are  Derman  [ 1 96 1 } , Roeloffs 
[1963,  1967],  Barlow,  Hunter  and  Proschan  [1963]  , McCall  [1969],  and  Zelen 
and  Feinleib  r 1969] . In  three  later  papers,  Keller  [1974],  Kirch  and  Klein 
[1974]  and  Pierskalla  and  Voelker  [1976],  models  are  presented  which  assume 
perfect  test  reliability.  In  these  papers,  optimal  mass  screening  sched- 
ules are  presented.  Keller  used  the  calculus  of  variations  to  characterize 
the  schedule.  Kirch  and  Klein  give  an  inspection  schedule  which  will  mini- 
mize expected  detection  delay  subject  to  a constraint  on  the  expected  number 
of  examinations  an  individual  would  incur  over  a lifetime.  Pierskalla  and 
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Voelker,  as  mentioned  before,  minimize  the  expected  disutility  per  unit 

time  for  the  population. 

•/ 

The  most  closely  related  work  to  the  results  in  this  paper  are  the 
works  of  Lincoln  and  Weiss  [1964]  and  Prorock  [ 1973 , 1976].  Lincoln  and 
Weiss  studied  the  statistical  characteristics  of  detection  delay  under 
the  assumption  that  the  times  of  examinations  form  a renewal  process  and 
that  the  probability  of  detecting  the  defect,  p(t),  is  a function  of  the 
defect's  age,  t.  They  derive  open  form  equations,  similar  to  renewal  type 
equations,  which  relate  the  density  functions  for  the  following  entities: 
the  probability  of  detection  at  a test  application  (p(t)),  the  time  until 
the  defect  becomes  potentially  detectable  and  from  this  time  the  forward 
recurrence  time  to  the  first  test,  the  probability  of  the  event  that  at  a 
particular  time  a test  occurs  and  all  prior  tests  hr.d  failed  to  detect 
the  defect,  and  the  detection  delay.  For  the  two  special  cases  where 
p(t)  is  a constant  and  where  p(t)  is  exponential,  the  moments  for  the  de- 
tection delay  are  derived  in  closed  form.  For  uniform  testing  intervals 
(and  general  p(*)),  the  distribution  and  moments  of  the  detection  delay 
are  computed.  They  also  studied  the  case  of  perfect  test  reliability. 

Prorok  has  established  properties  for  the  lead  time  in  the  case  of 
n successive  screening  tests  conducted  at  uniform  intervals.  He  also  de- 
rives the  proportion  of  preclinical  cases  detected,  both  at  a given  test 
and  over  all  n tests.  The  test  is  assumed  to  have  a constant  reliability 
not  necessarily  equal  to  one.  In  addition,  he  derives  an  estimator  for 
the  mean  lead  time.  This  estimator  utilizes  a mix  of  both  theoretical 
quantities  and  quantities  which  would  be  empirically  available.  In  the 
case  of  a perfectly  reliable  test,  his  estimator  requires  only  the  latter 
kind  of  input. 
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In  this  paper,  as  in  the  Pierskalla  and  Voelker  [ 1 976]  paper,  the 
model  does  not  seek  to  minimize  the  detection  or  inspection  delay  but 
rather  the  objective  function  is  (initially)  an  arbitrary  increasing  dis- 
utility function.  The  reason  for  choosing  a general  function  is  that  the 
disutility  experienced  upon  the  delayed  detection  of  a defect  may  well 
vary  in  a highly  nonlinear  way  with  the  length  of  the  delay. 

The  disutility  associated  with  a particular  detection  will  depend 
only  on  the  elapsed  time  since  the  defect's  incidence.  The  notation  D(t) 
will  express  the  disutility  incurred  if  detection  occurs  t units  of  time 
after  incidence.  D(* ) is  assumed  throughout  this  paper  to  be  a nonnega- 
tive increasing  function. 

In  those  cases  where  sup  D(s)  < °°,  there  is  a natural  relation  be- 

siO 

tween  the  notion  of  a utility  function  and  D(*);  namely, 

U(t)  = sup  D(s)  - D(t)  . 
s 2 0 

U(t)  is  a decreasing  function  expressing  the  disutility  avoided  by  de- 
tecting a defect  t units  of  time  after  its  incidence. 

In  the  next  section,  the  general  model  is  presented.  Section  2 con- 
tains the  analysis  of  the  case  when  the  probability  of  detection  is  a 
constant  over  all  values  of  elapsed  time  since  incidence.  The  long-run 
expected  disutility  per  unit  time  is  derived;  its  differential  qualities 
(with  respect  to  variations  in  the  test  reliability  parameter  and  testing 
frequency)  are  exhibited;  and  explicit  solutions  are  provided  for  various 
special  forms  of  D(* ) . Also  for  special  forms  of  D(*),  decision  rules 
are  presented  to  select  between  two  different  kinds  of  tests  which  differ 
with  respect  to  their  reliabilities  and  cost  per  application. 
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In  section  3,  it  is  assumed  that  the  test  will  detect  the  defect  if 
and  only  if  the  elapsed  time  since  incidence  exceeds  (or  equals)  a criti- 
cal threshold  T which  characterizes  the  tegt.  An  expression  for  the  long- 
run  expected  disutility  per  unit  time  is  derived,  and  decision  rules  are 
developed  to  select  between  two  alternative  test  types  which  differ  with 
respect  to  their  critical  threshold  T and  their  cost  per  application. 

Some  examples  for  linear  and  quadratic  disutility  are  given  in  the  final 
section. 

The  proofs  of  all  results  are  given  in  the  Appendix. 


1 . A General  Model  for  Time-Based  Mass  Screening 

This  model  was  developed  in  Pierskalla  and  Voelker  [1976]  and  a rela- 
tively simple  expression  for  the  objective  function  was  established.  The 


model  is  now  briefly  presented  (full  details  may  be  found  in  the  above 
paper) . 

It  is  assumed  that  the  jth  subpopulation  is  of  a fixed  size,  , and 
that  the  defects  arrive  according  to  a stationary  Poisson  process  with 
rate  NX  (N^  could  be  and  usually  is  a large  but  finite  number).  It  may 
be  somewhat  more  realistic  to  set  the  defect  arrival  rate  proportional  to 
the  number  of  defect-free  units,  rather  than  to  the  total  number  of  units 
in  the  population.  However,  it  is  also  assumed  that  no  defect  can  remain 
undetected  longer  than  T*  after  its  arrival,  even  without  any  screening 

“k 

tests  being  given.  Hence,  T is  an  upper  bound  on  the  expected  num- 

ber of  undetected  defects  in  the  jth  subpopulation.  If  it  is  the  case 
that  once  a defect  is  detected,  the  afflicted  unit  is  replaced  in  the 
population  with  a healthy  unit,  then  represents  a bound  on  the 
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expected  difference  between  the  number  of  healthy  units  and  the  total 
number  of  units  in  the  subpopulation.  For  small,  as  would  be  the 
case  for  a relatively  infrequently  occurring  disorder,  this  should  repre- 
sent no  difficulty.  In  this  paper,  therefore,  X ^ is  considered  to  be 
small,  which  is  consistent  with  the  examples  of  potential  applications 
which  are  mentioned. 

The  times  of  incidence  for  the  defect  in  the  jth  subpopulation  are 
designated  by  the  sequence  { S ^ } for  k =*  1,  2,...  The  screening  test  is 
assumed  to  be  administered  to  the  jth  subpopulation  at  the  times  1/r^, 
2/r^,  3/Tj,...  The  testing  frequency  r^  is  a control  variable.  The 
function  p^(» ) describes  the  reliability  characteristics  of  the  type  l 
screening  test,  i.e.,  how  likely  such  a test  is  to  detect  a defect  as  a 
function  of  the  defect's  age.  Note  that  p^(t)  = 0 for  t < 0. 

The  random  variable  S denotes  the  time  at  which  the  kth  defect 

r,i 

is  detected.  S depends  on  the  arrival  time  of  the  defect  (S  ) , the 
r ,X  j 

type  of  test  used  (A),  and  the  testing  frequency  (r^). 

Given  the  application  of  test  type  J l at  the  times  [l/r^,  2/r^,...}, 
the  disutility  incurred  by  the  kth  defect  is  where  D(*  ) is 

a strictly  increasing  function.  The  total  disutility  incurred  due  to 
those  defects  which  occurred  in  the  interval  [i/r^,  (i  + D/r^),  is 


<s?>' 


where  l^^ft)  is  an  indicator  function  defined  by 


1i/r(t)  ” 1[i/r,(i+l)/r](t) 


1 if  i/rs  ts  (i  + 1 )/r 


0 otherwise 
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The  following  theorem,  which  was  proved  in  the  previous  paper,  pro- 
vides an  expression  for  E[B  ] in  terms  only  of  the  disutility  due  to 

r , Z , i 

detection  delay  (as  -expressed  by  D(*))  and  of  the  reliability  of  test  type 

l (as  expressed  by  p (•)).  In  addition,  the  theorem  shows  that  E[B  .]  = 

l r , Z , 1 

E[Br,Z,0]  for  1 “ °>  l>  2’“- 


Theorem  A 


n 

00  r* 

f*  L i 


n-1 


E[Bi  i i1  = NX  S J J d<u)p  (u)  n [i-p  (u--^-)]du 

’ ’ n=l  , m=l 

n-1  J 


r . 
J 


for  i = 0,  1 , 2 , . . . 


By  defining  B^  as  the  long-run  expected  disutility  per  unit  time,  then 
r,Z 

by  Theorem  A 

r.  m-1 

S e[b->  , J 


r,Z  m . „ 

m-*“  1=0 


r,Z,iJ 


riE[Br  l 0]  = riNiXi  ^ I jD(u)p,(u)  n1  [1-p  (u-f-)Jdu 
J r,Z,0  J j j n=1  J i m=i  l r 

n-1  J 


r, 

J 


(The  factor  r^  enters  the  definition  to  convert  disutility  per  unit  testing- 
interval  into  per  unit  time.) 

Now  it  is  possible  to  formulate  a mathematical  program  for  the  problem 
of  selecting  testing  frequencies  and  test-types  for  each  of  Q different 
susceptibility  classes  which  together  comprise  the  whole  population.  These 
classes  may  differ  from  one  another  in  the  number  of  units  they  contain, 
in  their  defect  incidence  intensity,  and/or  in  the  cost  per  test  application 
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to  an  individual  for  a particular  type  of  test.  The  subpopulations,  how- 
ever, are  assumed  to  snare  a common  D(» ) function.  Also  it  should  be  noted 
that  the  test  used  may  be  a function  of  the  jth  subpopulation,  i.e., 

1 = i(J). 

The  multi-subpopulation  mass  screening  problem  subject  to  a budget 
constraint  is: 


(1) 

Minimize 

( r ^ , r 2 , . . • , r^ , ( 1 ) , • • 

t .... 

,,i(Q))  j-1 

such  that 

Q 

(2) 

S N.c.  „ , . . r,  £ b 
j=l  J J.Mj)  J 

(3) 

t—l 

It 

*•— » 

o 

A 

,r~ i 

u 

. Q 

(4) 

i(J>  €=£  J-1,... 

, Q 

where 


j 

b 


= cost  per  application  of  a test  of  type  £(j)  to  an 
individual  of  subpopulation  j , 

= number  of  units  (or  individuals)  in  subpopulation  j, 

= budget  per  unit  time. 


In  order  to  make  this  mathematical  program  even  more  comprehensive, 
it  is  possible  to  add  constraints  on  the  amount  of  testing  labor  available 
and  on  the  capacity  of  the  testing  facilities  in  terms  of  the  number  of 
arrivals,  the  frequency  of  testing  and  the  type  of  tests  used.  In  addi- 
tion, the  cost  of  false  positives  can  be  included  as  a part  of  the  test 
costs  in  inequality  (2). 
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For  example,  a constraint  on  the  total  labor  available  is: 


(5) 


1-1 


N 6 , . , . . r . s l, 
j j.-Kj)  J i 


and  on  the  total  testing  facilities  available  is: 


(6) 


£ N 

1-1 


. f . . , . r . ^ F, 
J J,4(j)  J l 


for  each  type  of  test  £(•)  used  over  the  subpopulations  being  tested, 
where 


6 ^(j)  = amount  °f  labor  needed  to  administer  test  type  2(j) 

to  an  individual  of  subpopulation  j , 

f.  = amount  of  testing  facility  time  needed  to  administer 

test  type  £(j)  to  an  individual  of  subpopulation  j, 

= total  amount  of  labor  available  to  administer  test 
type  £(• ) per  unit  time, 

F^  = total  amount  of  facility  time  available  to  administer 
test  type  £(• ) per  unit  time. 


This  mathematical  program  may  be  solved  by  any  of  several  nonlinear 
programming  methods.  However,  in  certain  special,  quite  important,  cases, 
more  explicit  formulations  of  the  objective  function  are  possible.  The 
previous  paper  by  Pierskalla  and  Voelker  [1976]  extensively  considered 
the  case  when  p^(t)  = 1 for  all  t > 0.  The  following  sections  are  con- 
cerned with  the  cases  of  constant  test  reliability  and  with  threshold 


test  reliability. 
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2.  Constant  Test  Reliability 

The  mass  screening  model  yields  interesting  results  for  a test  which 
has  a fixed  probability  p of  detecting  the  disorder  if  it  is  present  in 
an  individual.  Such  a model  would  arise  if  the  unreliability  of  the  test 
is  intrinsic  to  the  test  procedure  rather  than  partially  dependent  upon 
the  state  or  age  of  the  defect.  An  example  of  this  is  the  administration 
of  a Mantoux  test  for  tuberculosis  in,  say,  a population  of  grade  school 
children.  The  test  has  a small  but  relatively  constant  level  of  false 
negatives.  There  are  other  medical  tests  with  similar  characteristics 
(such  as  electrocardial  analysis). 

To  see  how  another  type  of  situation  with  constant  test  reliability 
could  arise,  consider  a production  process  which  is  subject  to  a randomly 
occurring  defect  which  degrades  its  performance.  Once  the  defect  occurs, 
the  level  of  degradation  of  the  process  remains  constant,  until  the  de- 
fect is  discovered.  Suppose  the  defect  is  such  that  each  item  produced 
has  probability  6 of  being  defective,  and  that  the  system  without  the  de- 
fect never  produces  defective  items. 

The  only  way  to  discover  the  existence  of  the  defect  in  the  produc- 
tion system  is  to  examine  an  item  produced  which  is  itself  defective. 

Now,  if  the  examination  of  an  item  is  expensive  (e.g.,  the  item  is  de- 
stroyed as  a result  of  the  inspection)  and  if  the  capacity  to  examine  a 
sequence  of  items  involves  a setup  cost  (say  a),  the  following  strategy 
might  uff  called  for:  At  specified  times  1/r,  2/r,.,.,  set  up  the  capacity 

to  examine  a sequence  of  items  and  examine,  say,  l items  at  each  of  those 
times.  The  times  1/r,  2/r,...  are  then  the  times  of  testing  and  the  sam- 
ple size  l € £ • (l,  2,  3,...}  specifies  the  test  type. 
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Assume  that  If  a defective  item  is  examined,  the  defect  is  always 
observed,  and  the  production  process  is  thereby  discovered  to  be  in  the 
degraded  state.  Hence,  a degraded  state  of  the  production  process  will 
go  undetected  at  the  testing  occasion  k/r  if  and  only  if  each  of  the  2 
items  sampled  at  time  k/r  is,  by  chance,  not  defective.  But  the  proba- 
bility  of  that  event  is  (1-6)  . Note  that  the  elapsed  time  t between 
the  entry  of  the  production  process  into  the  degraded  state  and  the  test 
time  k/r  does  not  affect  this  probability.  Hence, 

P4  = p2(t)  " 1 * (1  ‘ 6 

which  represents  the  probability  that  a test  (the  inspection  of  2 items) 
will  detect  a degraded  production  process. 

If  the  inspection  of  a single  item  has  cost  c,  then  to  set  up  and 
test  a sample  of  2 items  has  cost  c^  * a + 2c.  To  do  this  with  frequency 
r,  a budget  of  at  least  rc^  would  be  required. 

The  decision  maker  decides  how  large  a sample  of  parts  to  test  on 
each  testing  occasion,  and  how  frequently  to  schedule  those  occasions. 

If  there  are  a large  number  of  production  processes  involved,  they  may  be 
partitioned  according  to  their  likelihood  of  becoming  defective-- thereby 
forming  susceptibility  classes.  The  decision  as  to  sample  size  and  fre- 
quency would  then  have  to  be  made  for  each  class  so  as  to  stay  within  an 
t overall  budget  constraint. 

This  latter  problem,  several  susceptibility  classes  in  a population 
of  production  processes,  could  then  be  formulated  as  the  program  at  the 
end  of  the  previous  section.  For  present  purposes,  the  feature  of  inter- 
est in  this  example  is  that  p^(t)  does  not  depend  on  t. 
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Returning  to  the  more  general  development,  recall  that  . repre- 

r 

sented  the  long-run  expected  disutility  per  unit  time  if  test  type  i is 

applied  with  frequency  r . Exploiting  the  fact  that  p.(t)  is  a constant 

J *• 

(i.e.,  p^(t)  **  f°r  t 2 0)  and,  for  ease  of  analysis,  temporarily 
suppressing  the  subscript  j , 


rH-1 

— “ p 

(7)  Br  t - NXrp^  2 (1  - p^)n  J r D(u)du  . 

* * n«0 

n 

r 

In  order  to  avoid  confusion  with  the  general  case  of  B . given  in  the 

r ,* 

preceding  section,  it  is  convenient  to  use  a new  notation  for  (7).  Conse- 
quently, let  p be  denoted  by  p and  let 

AJ 

n+1 

00  * 

(8)  C(r,  p)  » N\rp  72  (1  - p)n  r D(u)du  . 

n-0  J 

n 

r 

Therefore,  (7)  and  (8)  represent  the  long-run  expected  disutility  per 
unit  time  if  a test,  of  a type  which  is  characterized  by  a reliability 
which  does  not  vary  with  time  since  the  defect's  incidence,  is  adminis- 
tered with  frequency  r. 

A question  which  is  examined  next  is  -how  do  changes  in  r and  p 
affect  C(r,  p).  After  that  explicit  solutions  are  given  for  C(r,  p) 
when  D(* ) takes  certain  simple  forms.  And  lastly,  some  general  rules  are 
indicated  for  selecting  between  a particular  kind  of  test  and  a more  ex- 
pensive but  more  reliable  alternative  test  when  D(* ) takes  certain  forms. 

The  following  lemma  indicates  that  C(r,  p)  behaves  consistently  with 
one's  expectations  as  r or  p varies. 
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Lemma  1:  If  D(* ) Is  a strictly  increasing  function, 

(9)  o , and 

dp 

(10)  ac-fr  p)  < 0 . 

o r 


Note  that  for  D(* ) nondecreasing,  the  inequalities  in  (9)  and  (10) 
still  hold  but  not  strictly. 

From  this  lemma,  as  anticipated,  when  p increases,  the  expected  dis- 
utility decreases.  Similarly  as  r increases,  the  interval  1/r  between 
tests  decreases  and  the  expected  disutility  decreases.  Consequently,  as 
better  test  types  are  used  or  the  tests  are  more  frequently  applied,  the 
value  of  such  changes  in  terms  of  reduced  disutility  versus  the  costs  of 
the  changes  can  be  assessed  and  the  tradeoffs  evaluated. 

It  is  easy  to  compute  the  Hessian  for  C(r,  p)  when  D(* ) is  differ- 
entiable. 


^-N\pS  qV3  f (n+l)2D'(— ) - n2D'  <S)“| 

dr  n-0  L r r J 

2 CO 

— j » NXr  Tj  n(n-l)qn  [W(n,  r)  - W(n  - 1,  r)] 
dp  n-2 

2 ® 

m NX  Ti  [l  - (n  + l)p]qn  1[W(n,  r)  - A(n+  1,  r)  + A(n,  r)] 
3pdr  n-0 

rH-1 

where  q - 1 - p;  W(n,  r)  - (*  r D(s)ds;  A(n,  r)  - — D(— ) . 

J n r r 
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& C 9 C 

Note  that  if  D' (• ) is  increasing,  then  — t 2 0 and  — j 2 0.  Hence, 

Sr  Sp 

along  coordinate  directions  C(r,  p)  is  convex. 

Simple  expressions  for  C(r*  p)  can  be  given  when  D(* ) is  specialized 
to  a geometric  and  an  exponential  function,  respectively.  Since  these  two 
types  of  functions  are  reasonably  general  and  are  intuitively  appealing 
for  disutility  functions  (in  the  sense  of  an  increasing  rate  of  disutility 
the  longer  the  defect  is  undetected),  they  can  be  quite  useful  as  realistic 
approximations  in  applications . 


Lemma  2:  If  D(t) 


atm  for  a > 0 and  m 2 

00 


C(r,  p) 


aN  \p  y> 
(m+  l)rm  n*l 


1,  then 

ttrt-1  n- 
n pq 


1 


In  the  case  that  m is  £ positive  integer,  then  the  summation  in 
Lemma  2 above  is  just  the  moment  generating  function  of  a geometric  ran- 
dom variable  with  parameter  p and  C(r,  p)  can  be  expressed  by 


(11) 

where 

For  m 


C(r,  p) 


aN\pl(mfl)(0  : p) 
(m  + l)rm 


i|r(t  : p)  » — r and  ^m^(0  : p)  = 

1 - qeC  dt"1 


t = 0 


1,  C(r,  p)  - (~~^) , and  for  m - 2,  C(r,  p)  - — £■  \ l + 6 : 2&P  1 

p 3r  — p 


Lemma  3:  If  D(t)  ■ geat  for  a,  g > 0,  then 

C(r,  p)  - (3N\pr  (e3/,r  - l)/a(l  - qe3^) 


for  r > 


- a 

log(q) 
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Lemmas  2 and  3 can  provide  a means  to  select  between  two  alternative 
kinds  of  tests  which  differ  with  respect  to  reliability  of  detection  and 
cost  per  application.  Let  test  #1  have  cost  per  application  c^  and  relia- 
bility p^.  The  corresponding  parameters  for  test  #2  are  c^  and  p2>  If 
test  #1  is  administered  with  frequency  rc^/c^,  and  test  #2  administered 
with  frequency  r,  both  testing  regimes  will  consume  equal  quantities  of 
the  budgeted  resource;  viz,  Nrc2>  If  C(rc2/c^,  p^)  £ C(r,  p^)  for  all 
r 2 0,  then  the  expected  disutility  per  unit  time  will  be  less  with  test 
#1  at  all  levels  of  budget  Nrc0.  That  is,  if  test  #2  is  being  used  with 
frequency  r,  the  expected  disutility  can  be  decreased  without  any  addi- 
tional allotment  of  budget,  simply  by  switching  to  test  #1  and  testing 
as  frequently  as  the  budget  permits. 

Suppose,  for  example,  that  D(t)  = atm  for  m a positive  integer.  Then 


C(rc2/c1,  pL)  £ C(r,  p2> 


if  and  only  if 

(12) 


(c1/c2)m  £ 


P2^  (°  : P2> 

P^T(oTp~7 


If  the  inequality  in  (12)  is  reversed,  then  test  #2  is  preferred  to  test 
#1 . Notice  that  this  rule  does  not  depend  on  r,  and  hence  does  not  depend 
on  the  budget  contribution  Nrc,,. 
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3 . Threshold  Test  Reliability 

In  the  previous  section  the  reliability  of  the  test  depended  only  on 
factors  intrinsic  to  the  test  itself  and  did  not  depend  at  all  on  the 
elapsed  time  since  incidence  at  the  time  of  the  test.  In  this  section  a 
special  form  for  p(t)  is  considered  which  is  very  different  from  the  case 
of  constant  test  reliability.  Here  the  test  reliability  is  zero,  if  the 
elapsed  time  since  the  defect's  incidence  is  less  than  T,  and  that  at 
elapsed  time  T,  the  reliability  jumps  to  one.  That  is,  p(t)  « 1^  ^(t) 
where  the  number  T is  a characteristic  of  the  type  of  test  chosen. 

The  primary  results  in  this  section  are  a simple  characterization  of 
B . when  p (•  ) * lr_  .,(•),  and  rules  which,  in  some  cases,  will  permit 

r>*  l l » / 

selection  between  two  tests  which  differ  in  their  reliability  (i.e.,  in 
their  detection  threshold  T)  and  in  their  cost  per  application. 

Since  the  number  T^  completely  specifies  the  function  p^(t)  “ 1^T  ^(t), 

the  long-run  expected  disutility  per  unit  time  is  expressed  in  terms  of 
r and  T rather  than  in  terms  of  r and  the  test  type  index  i.  Suppressing 
the  subpopwlation  index  j and  the  test  index  l,  define 


(13) 

When  T 
of  PA(t)  - 1 


A(r * T>  3 Br,4  When  pA(,)  " 1[Tj4>-)(,)  • 

0,  this  definition  is  equivalent  to  perfect  test  reliability 
for  t > 0 and 

l 

A(r , 0)  * rNX  f r D(u)du  . 

v 

0 


The  next  theorem  evaluates  A(r,  T)  for  all  values  of  T i 0.  Of 
course,  the  "blind  period"  of  the  test  from  0 to  T does  not  delay  the 
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detection  of  each  arriving  defect  exactly  T units  of  time.  The  amount  of 
delay  depends  on  the  interplay  among  the  time  of  arrival  of  the  defect, 
the  testing  schedule  {l/r,  2/r,...},  and  the  magnitude  of  T. 

Theorem  1:  If  p(*  ) = 1^  „)(*  ) for  some  T a 0,  then 

\ 

(14)  A(r,  T)  = NXr  j r D(T+u)du  . 

0 

Suppose  the  decision  maker  has  two  kinds  of  tests  available,  and  he 
must  choose  one  of  them  for  implementation  in  a mass  screening  program. 
Suppose  the  first  kind  of  test--call  it  test  #l--has  sensitivity  charac- 
terized by  the  "time-until-detectability"  threshold  T^ . That  is,  this 
test  will  detect  a disorder  if  and  only  if  the  disorder  has  been  present 
for  a length  of  time  T^ . Let  c^  > 0 be  the  cost  per  application  (to  an 
individual  unit)  of  this  kind  of  test.  For  the  second  kind  of  test  under 
consideration,  test  #2,  let  T^  and  c^  be  the  corresponding  parameters. 

Assume  test  #1  is  better  in  the  sense  that  T^  < T£.  If  c^  £ c^, 
there  is  no  decision  to  make;  test  #1  should  be  used.  Therefore,  assume 

c I3"  c2' 

If  the  exact  shape  of  the  function  D(* ) is  known,  Theorem  1 can  be 
used  to  decide  which  test  to  use  for  each  possible  level  of  budget.  Let 
b be  the  budget  per  individual  in  the  (homogeneous)  population.  Then  the 
use  of  test  #1  will  permit  a testing  frequency  of  b/c^,  and  the  use  of 
test  #2  permits  frequency  b/c2*  To  decide  which  test  to  use,  compare  the 
expected  disutilities  per  unit  time  assuming  a fully  allocated  budget, 
i.e.,  compare  A(b/c^,  T^)  and  A(b/c2>  T2).  With  D(* ) known,  these  quanti- 
ties can  be  evaluated  explicitly  by  Theorem  1 and  compared. 
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Note  that  it  was  assumed  that  the  entire  budget  should  be  allocated. 
This  assumption  was  unnecessary  because  when  D(» ) is  an  increasing  function 

l 

^ A(r , p)  - NX  \ j rD<T  + u)  - D(T  + ^)du]  < 0 . 

“ 0 


Therefore,  the  greater  the  testing  frequency,  the  lower  the  expected  dis- 
utility per  unit  time.  For  this  reason,  the  budget  should  be  entirely  al- 
located. 

When  the  exact  form  of  the  disutility  function  is  not  known,  Theorem 
1 does  not  directly  differentiate  between  test  #1  and  test  # 2.  However, 
the  two  following  theorems  will  permit  such  a determination,  at  least  for 
certain  relative  configurations  of  the  budget  (per  unit  population),  and 
the  relative  sensitivity  and  cost  differential  c2  - of  the  tests. 

Specifically,  Theorem  2 will  show  that  for  any  (increasing)  disutil- 
ity function  test  #1  is  indicated,  if  the  budget  (per  unit  population) 


exceeds  (c^  - c^)/ (T^  - T^) . On  the  other  hand.  Theorem  3 shows  that  for 
a convex  increasing  disutility  function,  test  #2  is  better  if  the  budget 


is  less  than 


c - c 
1 2 

2(T2  - Tl)  * 


A decision  rule  for  the  case  where  the  budget  falls  between 


Cl-C2 


■~j= has  not  been  found  for  general  disutility  functions. 

' 2 * V 


Just  how  the  statements  of  Theorems  2 and  3 are  translated  into  the 


above  decision  rules  is  explained  after  the  statements  of  the  respective 


theorems . 
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Iheorem  2:  Given  D(* ) a strictly  increasing  function,  and 

c^  > c then  - T > (c^-c^)/rc^  provides  a sufficient  condition  for 


(15) 


A(rc2/Cl,  Tl)  < A(r , T2> 


and  hence  for  test  #1  to  be  preferable  at  the  per  unit  population  budget 
level  of  rc2* 

The  hypothesis  of  Theorem  2 is  equivalent  to 


(16) 


r > 


CfC2 

c2(T2  ' Tl) 


If  test  #1  is  used  with  frequency  (as  on  the  left  side  of 

(15)),  the  budget  required  would  be  rc2  which  is  the  same  budget  required 
to  execute  test  #2  with  frequency  r (as  on  the  right  side  of  (15)). 

Hence,  translating  (16)-  into  its  equivalent  form  in  terms  of  budget 
(per  unit  population),  designated  by  b, 


b = rc2  > (c1  - c2)/(T2  - Tx)  . 

With  its  hypothesis  in  this  form,  the  theorem  provides  a lower  bound  on 
the  budget  which  is  a sufficient  condition  for  test  #1  to  entail  lower 
expected  disutility  per  unit  time  vis-a-vis  test  #2,  were  the  two  tests 
scheduled  at  their  respective  maximal  frequencies  rc^/c^  and  r consistent 
with  the  budget  rc2*  The  following  lemma  is  needed  for  the  proof  of 
Theorem  3 and  is  recorded  here  for  general  interest. 


Lemma  4:  If  f is  a convex  function,  then 


1 Tt+y  1 

“ j f(s)ds  * y[f(t)  + f(t  + y)]  . 

t 


(17) 


-23- 


Figure  1 

This  Figure  illustrates  Theorems  2 and  3.  Each  point  on  the 
sloping  line  corresponds  to  a level  of  budget.  Regions  where 
each  test  is  preferable  are  indicated. 
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Iheorem  3:  If  D(* ) is  convex  and  increasing,  c^  > c2>  T^  < , and 

T2  - TL  s (c^^  - c2)/2rc1 , then 

(18)  A(Clr/c2,  T2)  s A(r,  . 


The  hypothesis  of  Theorem  3 is  equivalent  to 
(19)  r £ 


C1  ' C2 


2(T2  - T1)c1 


If  test  #2  is  used  with  frequency  c^r/c2  (as  on  the  left  side  of 
(18)),  a budget  (per  unit  population)  of  c^r  would  be  required.  The 
same  budget  is  required  to  use  test  #1  with  frequency  r (as  on  the  right 
side  of  (18) ) . 

Letting  b represent  the  budget  per  unit  population,  (19)  is  equiva- 
lent to  the  following  constraint  on  the  budget: 

C1‘C2 

b Clr  S 2(T2  - Tx)  • 


Theorem  3 then  indicates  the  superiority  of  test  #2  at  such  levels 
of  budget. 


4.  Some  Examples  of  Disutility  Functions 

Suppose  a production  process  is  subject  to  a randomly  occurring  de- 
fect. Although  production  appears  to  proceed  normally  after  the  incidence 
of  the  defect,  the  product  produced  is  thereafter  defective  to  an  extent 
which  remains  constant  until  the  production  process  is  returned  to  its 
proper  mode  of  operation.  The  only  way  to  learn  if  the  production  process 
is  in  this  degraded  state  is  to  perform  a costly  test.  Now,  if  a test  de- 
tects the  existence  of  the  degraded  mode  of  production  t units  of  time 
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after  its  incidence,  the  harm  done  will  be  proportional  to  the  amount  of 
defective  product  (unknowingly)  produced  which,  in  turn,  is  proportional 
to  t.  Hence  D(t)  * at  for  some  a>  0. 

Another  example  where  a linear  D(* ) function  may  be  appropriate  would 
be  for  the  periodic  inspection  of  an  inactive  device  (such  as  a missile) 
stored  for  possible  use  in  an  emergency.  If  t is  the  time  between  the  in- 
cidence of  the  disorder  and  its  detection,  the  disutility  incurred  may  be 
proportional  to  the  probability  that  the  device  would  be  needed  in  that 
time  interval.  If  such  "emergencies"  arise  according  to  a Poisson  process 

with  rate  i*  , then  the  probability  of  an  emergency  in  a time  interval  of 
t 

length  t is  1 - e , which  for  small  is  approximately  ^ t.  Hence,  if  b 
is  the  cost  incurred  should  there  be  an  emergency  while  the  device  is  de- 
fective, and  if  i*  is  the  (small)  arrival  rate  of  emergencies,  then 
D(t)  * bn  t . 

A quadratic  disutility  could  arise  in  the  following  situation.  Sup- 
pose the  magnitude  of  a randomly  occurring  defect  increases  linearly  with 
time  since  the  occurrence  of  the  defect.  For  example,  the  magnitude  of 
the  defect  might  be  the  size  of  a small  leak  in  a storage  container  for 
a fluid,  and  as  fluid  escapes  the  leak  gets  larger.  Further  suppose  that 
the  harm  done  accumulates  at  a rate  proportional  to  the  magnitude  of  the 
defect.  Hence,  the  quantity  of  fluid  lost  (at  least  initially)  increases 
the  longer  the  defect  exists,  and  the  rate  of  fluid  loss  is  proportional 
to  the  size  of  the  leak. 

Let  the  size  of  the  leak  (as  measured  by  rate  of  fluid  loss),  at 
time  s since  the  leak's  incidence,  be  cs . Then,  if  the  defect  is  detected 
at  time  t since  incidence,  the  disutility  incurred  (fluid  lost)  is 

D(t)  - J csds  = 2Ct^  ’ 

0 
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Examples  of  nonlinear  disutility  functions  for  a few  animal  and  human 
diseases  are  already  known.  Although  data  may  not  be  available  currently 
to  support  the  utilization  of  particular  models,  that  should  not  deter  the 
development  of  such  models.  It  is  reasonable  to  expect  that  in  many  ap- 
plication areas  in  the  future  much  additional  data  will  become  available; 
for  example,  data  concerning  the  stochastic  pattern  of  a particular  dis- 
ease's development.  Furthermore,  these  models  can  serve  as  a guide  to 
the  kinds  of  data  which  should  be  gathered. 

It  is  also  reasonable  to  expect  the  future  development  of  improved 
testing  technologies  capable  of  detecting  a disorder  in  a much  earlier 
stage  of  development  than  is  now  the  case.  Such  innovations  will  make 
screening  programs  more  attractive  and,  consequently,  make  more  impor- 
tant the  analytical  tools  to  design  such  programs  intelligently. 
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Proof  of 


where 


Now  let 


Then,  W 

n 


(Al) 


Appendix 


Lemma  1:  Letting  q - 1 - P, 


A 


C(n,  p) 


NXr  £ q'VCn)  - NXrp  £ 
n-0  n-0 


nq11  Hf(n) 


n+1 

W(n)  - J r D(s )ds  . 

n 

r 


vo-wo 

Vn  “ Wn'Vl  ’ n “ l* 


J-0 


and 


i n 

n-1  r' 


C(r,  p)  - NXr  £ q”  £ V - N\rp  £ nq  ‘ £ V 

OP  n j _A  J A 1_A  J 


n-0  j *0 


n**0 


j-0 


NXr 


Z if  S ,‘.,I  nq”'1] 

j-0  J w n-j  n-j 

CO 

- -NXr  £ jqj"V 


j-0 


. i J±L 

- NXr  £ jqJ’1  i r D(s)ds  - J r D(s)dsJ  , 
j_1  hi  1 


£ qn 

n-j 


q - p Li  nq 
n-j 


n-1 


•jq 


J-i 


since 


qJ"l[q/p  - j - q/pl 
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Because  D(» ) Is  strictly  increasing,  the  right-hand  side  of  (Al)  is  nega- 
tive, and  C(r,  p)  < 0. 

OP 


dr 


C(r,  p) 


NXp  £ 
n«0 


nfl 

q11  J r D(s  )ds 

n 

r 


- N\p  2 
n-0 


nf  1 
r 


D ( 


n+  1 
r 


) - 


= NXp  £ qnT  | r D(s  )ds  " “ (d  (^-i)  - D <7>)  -7°(il7i)j  < 0 • 
n— 0 k_ 

n 

r 


This  inequality  follows  from  D(« ) increasing,  through  the  relations 
,n  + 1, 


and 


D(- 


r -)  - D(^)  > 0 


nfl 


nfl 


f r D(s )ds  < [ r D(^^)ds  -7D(^-^)  - 0 . 

j r r j r r r 


n 

r 


Q.E.D. 


Proof  of  Lemma  2: 


at I 

C(r,  p)  » aNXrp  £ q11  f r tmdt 

n-0  J 
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■ aN*e--  | s qV+i)"^1 

(m  f 1 ) r”1  “ n-0 


90 
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2^  q n 

n=0  J 
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n-1  mfl 
q n 
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q n 
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Q.E.D. 
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Proof  of  Lemma  3: 


\ r k r^ak/r,  a/r  . 

C(r,  p)  - NXp  2_i  q — e (e  - 1) 

k-0  3 


NXpr  , a/r 


.a/r  . x r\ 

(e  - 1)  L (qe  ) 


a/r.k 


00 

_ 3/  IT  Ic  s/  IT 

The  geometric  series  (qe  ) converges  if  and  only  if  qe  < 1, 

k-0 


i.e.,  if  and  only  if 


log  q 


Therefore,  for  r> 


log  q 


NXor  , a/r  . * 1 

C(r,  p)  - — f"(e  -1)  S77  • 

1 - qe 


Proof  of  Theorem  1:  The  following  conventions  are  used  below: 


II  x - 1 for  n < 1 and  [x]  is  the  largest  integer  not  exceeding  x. 
1-1  1 


From  the  definition  of  B 

r ,x 


00  (*~  n-1  r -1 

(A2)  A(r,  T)  - NXr  S r D(u)lf  . (u)  n | 1 - lr  , (u  - ^ ) 

n-1  J , 1 ’ ; m-1  *”  1 * ’ r J 


[ ru]  r 

«“>-  n,Ll-l[T.-)<u-r>J 


Note  that  for  u € [ (n  - l)/r,  n/r),  n-1  - [ru] 
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Therefore,  (A2)  becomes 


" r“ 

A(r,  T)  » NXr  Ti  J r D(u)1jt  ^ (u) A (u)du 
n“l  __  , ’ 


Now  it  will  be  shown  that 


A “ 1[0>T  + ~)  * 


Select  x 2 T + - . Then  [rx]  2 1 and  by  (A3) 


[ rxj  p , 

■w  ■ ”,  L1'l[T.-)<x"?)]  5 

m»I 


But  x-  — 2 T,  therefore  the  above  relation  shows 
r 

(A6)  A(x)  =»  0 for  x 2 T + -^  . 

Now  select  y€[0,  T + -^  ) . If  [ry]  = 0,  (A3)  and  the  convention  on  de- 
generate products  gives  A(y)  » 1 . On  the  other  hand,  if  [ry]  > 0,  then 
the  relation  y-~s  y*“<  T for  m “ 1>  2,...,  [ry]  is  not  vacuous  (i.e., 
the  set  of  allowable  indices  m is  not  empty)  and  neither  is  the  consequent 


relation 


1 [T  „ ) <y  - ~ > ■ 0 for  b - 1 , 2 , [ ry]  . 


And,  hence  by  (A3) 


A(y)  = 1 for  y €[0,  T + ±)  . 


(A5)  now  follows  from  (A6)  and  (A7).  Substituting  (A5)  into  (A4) , 
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f co  a 

A(r,  T)  - NXr  S j r D(u)l[x>oo)  («)1[0  T+  l}  (u)du 

n"*  n-1  r 

r 

- NXr  J D(u)l[T>T+l)(u)du 

0 * r 
l 

- NXr  J rD(T+  s)ds  . 

0 


Q.E.D. 


Proof  of  Theorem  2: 

Suppose  t“*D(t)  Is  the  underlying  disutility  function.  Then,  by 
Theorem  1,  (15)  is  equivalent  to  showing 


(A8)  J 2 D(T1  + s)ds  < r j rD^T1+ (^-T^+s^ds  . 

1 0 


Let 


D (t)  - D(TX  + t)  and  T - T2  - ^ , 


then  (A8)  becomes 


(A9) 


rc-  -c  /rc 

— J D (s)ds  < r J r D*(T  + s)ds  . 


Since  D(» ) is  an  arbitrary  increasing  function,  so  is  D*(»);  consequently, 

•ff 

the  asterisk  on  D (•  ) may  be  dropped. 

To  prove  the  theorem,  it  suffices  to  verify  (A9) . By  hypothesis, 

T > (cx  - c2)/rc2 

which  implies 

for  u € [0, 


T > u(ct  - c2)/c2 
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which  implies 

(A10)  T + u>  uc^/c2  for  u € [0,  . 

Now, 

/ I 1 

rc_  »c  /r c_  » p 

j 1 Z D(s)ds  - r j D(c1s/c2)ds  < r J D(T+s)ds  . 

1 0 0 0 

The  inequality  follows  from  (A10)  and  the  fact  that  D(* ) is  increasing. 

Q.E.D 


Proof  of  Lemma  4: 

Let 

h(t+s)  - f (t)  [f(t+  y)  - f (t)]  , s € [0,  y]  . 

Since  f is  convex, 

f(t+s)  £ h(t  + s)  , s € [ 0 , y]  . 

Hence , 

pt+y  l pt+y  py 

— J f(s)ds  s - J h(s)ds«—  j h(t  + s)ds 

t t 0 

y “ 

-“J  * f(t)  +“(  f(t  + y)  - f(t))]  ds 
0 

- f(t)  + \ [f (t  + y)  - f(t)](y2/2) 

y 

- -j  [f(t  + y)  + f ( t )]  . 


Q.E.D 
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Proof  of  Theorem  3: 

Theorem  1 renders  the  inequality  (18)  equivalent  to 


(All) 


c r pC./c  r 

— J D(T2  + s)ds  s r j rD(T1  + s)ds 

2 0 0 


As  in  the  proof  of  Theorem  3.6,  let 


D (t)  - D(T^ + t)  and  T - - T . 


(All)  then  becomes 


(A12) 


c r rc,/c  r - 

~ J D (T+  s)ds  s r J D*(s)ds  . 


Note  that  D(* ) being  an  arbitrary  convex  increasing  function  implies  D*(* ) 
is  an  arbitrary  function  of  the  same  class.  Hence,  the  asterisk  may  be 
dropped.  Also,  let  a - and  X “ r ‘ Then>  (A12)  can  be  written 


(A13) 


1 »ax  x 

-J  D(T+  s)ds  s ■“  J D(s)ds 


Using  Lemma  4 
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0 T 
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T 
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(A14) 


T 

- (1/2T)  J G(0)ds 
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where  G(s)  = D(T  - s)  + D(T+  ax+  s) . Now  D(*  ) convex  and  increasing  implies 
G(* ) is  convex  and  increasing  in  s.  Therefore,  when  s a 0,  G(0)  £ G(s) 
which,  along  with  (A14) , implies 

i <»ax  i i*  ^ i r ^ 

D(T+s)ds£-^J  G(0)ds  £ J G(s)ds 
0 0 0 

T 

(A15)  “ IT  I [D(T- s)  + D(T+ax+s)]ds  . 
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„2T+ax 
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This  last  inequality  follows  from 
T+ax  , T 

D(s)ds  s J [D(T  - s)  + D(T  + ax+s)]ds 
T 0 

which  follows  from  (A15).  By  hypothesis,  1 s x(c^  - c^)/2c^  which  implies 
(A17)  2T  + ax  s x . 


,2T+ax 


1 A A -l  (,  +-  X f OA  -*  A A 

— I D(s)ds  ■ — I D(s)ds  + — J D(s)ds 


2T+ax 
2 T+ax 


/ 2T  + ax\  1 I*  " 

D(s)ds 


(A18) 
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Hie  last  inequality  follows  from 

,2T+ax 
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2T+  ax 


D(s)ds  s x - 2T  - ax  J D<s> 


ds 


2T+ax 


which  follows  from  (A17)  since  D(* ) is  increasing.  The  proof  is  completed 
by  using,  in  turn,  (A16)  and  (A18)  to  obtain 

T+ax 


. pax  pT+ax 

-J  D(T+s)dS--l  D(a)ds 


<1 
T 

p2T+ax  P x 

s IFTS  J D<3)ds  s x J D<3)ds 


q.e.d. 
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20. 

In  both  cases,  the  defect  arrival  process  is  Poisson. 

In  the  first  case,  the  reliability  of  the  test,  i.e.,  the  probability 
that  the  test  will  detect  the  disease,  is  dependent  only  on  the  type  of  test 
being  administered  and  is  constant  over  all  values  of  elapsed  time  since  the 
incidence  of  the  defect.  The  long-run  expected  disutility  per  unit  time  is 
derived  and  variations  of  the  disutility  with  regard  to  test  reliability  and 
testing  frequency  are  presented.  In  addition,  explicit  solutions  are  provided 
for  various  special  forms  of  the  disutility  function. 

In  the  second  case,  the  test  will  detect  the  defect  if  and  only  if  the 
elapsed  time  since  incidence  exceeds  a critical  threshold  T which  character- 
izes the  test.  An  expression  for  the  long-run  expected  disutility  per  unit 
time  is  derived  and  decision  rules  to  select  between  two  alternative  test 
types  and  costs  are  presented. 
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